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Abstract 

We use the worldline path-integral approach to the Bern- 
Kosower formahsm for developing a new algorithm for calculation 
of the sum of all diagrams with one spinor loop and fixed num- 
bers of external and internal photons. The method is based on 
worldline supersymmetry, and on the construction of generalized 
worldline Green functions. The two-loop QED f3 - function is 
calculated as an example. 
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1. Introduction 

In 1992, Bern and Kosower [1] used string theory to derive a new for- 
malism for the calculation of one-loop amplitudes in ordinary quantum field 
theory which is equivalent to Feynman diagrams [2], but leads to a sig- 
nificant reduction of the number of terms to be computed in gauge theory 
calculations. This property was then successfully exploited to obtain both 
five-point gluon [3] and four-point graviton amplitudes [4]. 

Following this, Strassler [5] showed that, at least for scalar and spinor 
loops, the resulting integral representations can be derived in a more ele- 
mentary way. In this approach, one writes the one-loop effective action as a 
(super)particle path-integral, and evaluates this path-integral in analogy to 
the Polyakov path integral, i.e. using worldline Green functions appropriate 
to a one-dimensional field theory on the circle. 

This reformulation turned out to be useful for various calculations of one- 
loop effective actions [6, 7, 8] and yielded, in particular, a new method for the 
calculation of the inverse mass expansion which is nonrecursive, manifestly 
gauge invariant, and suitable to computerization [7, 9]. 

Progress has been made along different lines to generalize the Bern- 
Kosower formalism beyond one loop, using methods either based on the 
calculation of higher genus string amplitudes [10], on the use of a separate 
worldline path integral for every internal propagator [11], or on a stringlike 
reorganization of standard Feynman parameter integrals [12]. A Hamilto- 
nian approach has also been considered [13]. 

Recently, we have proposed [14] a multiloop generalization of Strassler's 
approach, based on the concept of worldline Green functions for multiloop 
diagrams. Those Green functions have been explicitly constructed for the 
general two-loop graph, and for a loop with an arbitrary number of propaga- 
tor insertions. If used with global proper-time variables, this allows to derive 
integral representations combining whole classes of Feynman diagrams into 
compact expressions. 

While knowledge of the worldline Green functions is, in principle, suffi- 
cient to treat arbitrary scalar diagrams, more work has to be done to obtain 
the final integral representations for multiloop amplitudes in general quan- 
tum field theories. 

In the present paper, we take up the study of quantum electrodynamics, 
and consider a class of amplitudes which is the simplest one for our purpose, 
namely the N-photon amplitude with a single spinor loop. 
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2. One-Loop Amplitudes 

First let us shortly review how one-loop calculations are done in this 
formalism [5, 7, 14]. For calculation of the one-loop effective action induced 
by a massive spinor-loop in a background gauge field, one would start with 
the following worldline path integral representation [15, 16, 17, 18, 19]: 



T[A] = -2r^e-"''^tT[VxVi; 
Jo T J 

xexp 



/■^ /I 1 • 

Jo ^4 Z 



Here the a:;^(T)'s are the periodic functions from the circle with circumfer- 
ence T into D - dimensional spacetime, and the 'ip^[Tys their antiperiodic 
Grassmannian supersymmetric partners. In the nonabelian case - which we 
will not consider in this paper - path ordering would be implied. 

For calculation of the effective action, one first splits the coordinate path 
integral into center of mass and relative coordinates. 



rdry^ir 
Jo 



J dxo J Vy 
x'o + y'ir) 
. 



One then Taylor-expands the external field at xq, and evaluates the path 
integrals over y and ip by Wick contractions, as in a one-dimensional field 
theory on the circle. The Green functions to be used are those adapted to 
the (anti-) periodicity conditions. 



{r{n)r{r2)) 



\g^''Gp{n,T2] 



-9' 



,fJ,U 



n - T2 



[n-r2 



-^^''sign(Ti - T2] 



(3) 



Our normalization is such that for the free path integrals 
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f /"-^ 1 

J T^ip exp^— J dT—iptp^ = 1 

(4) 

The result of this evaluation is the one-loop effective Lagrangian C{xq). 

One-loop scattering amplitudes are obtained by specializing to a back- 
ground consisting of a finite number of plane waves. This amounts to the 
same thing as defining integrated vertex operators 



f 

Jo 



dr [x^Cu — 2iip^ip'^ ki^e^]ex]i[ikx[T)] (5) 

/o 

for external photons of definite momentum and polarization, and calculating 
multiple insertions of those vertex operators into the free path integral. 

Using the worldline superfield formalism of [22, 18], this calculus may be 
cast into manifestly supersymmetric form. Eq.( 1) then becomes 



r[A] = -2 r^e-'^ /i)Xexp-^o^^^^^4-^^^^^-^^^^^(^)l, 
Jo T J 



(6) 



where 



D - ^ 9^ 



dee = 1 

(7) 

The photon vertex operator is rewritten as 

- r drdee.DX^'expiikX], (8) 
Jo 

and the worldline propagators may be combined into a superpropagator 

G{n,e,-T2,e2) = GB{ri,T2) + e,e2GF{n,T2). (9) 

From this superfield formalism, it is not difficult to derive the following im- 
portant "substitution rule" [1, 5, 6]: 
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Evaluation of the bosonic path integral in general leads to an expression 
consisting of an exponential factor exp Z)i<j GB{Ti,Tj)piPj^, multiplied by a 
polynomial in the first and second derivatives of Gb, 



Gb{ti,T2) = sign(Ti - T2) - 2 



Tl - T2 



(10) 



(here and in the following, a "dot" denotes differentiation with respect to 
the first variable). The G^'s can always be eliminated by partial integra- 
tions on the worldline, and once this has been done, all contributions from 
fermionic Wick contractions may be taken into account by replacing every 
closed cycle of Gs's appearing, say (^^(tji , t,2)G's(t,2, r^j) • • • G's(t,„,t,J, by 
its supersymmetrization: 



G's(t,j,t,JG's(t,2,t,3) • • • G's(t,„,t,J GB{n^,n^)GB{n^,n^) ■ ■ ■ GB{n^,n^) 

- Gf{t,^,t,^)Gf{t,^,t,^) ■ ■ ■ Gf{t,^,t,^) 

(11: 

Note that the result would vanish if supersymmetry was not broken by the 
boundary conditions. This substitution rule effectively replaces the calcula- 
tion of Dirac traces. 

3. Scalar Multiloop Amplitudes 

In our treatment of the two-loop scalar diagram [14], the starting point 
had been to insert into the free one-loop scalar path integral 



/ dTa / dn {4){x{Ta))4){x{n))) . 
Jo Jo 



a scalar propagator 



:i2i 



(13) 



This propagator was then written in the Schwinger proper-time representa- 
tion. 



{<f>{x{r^))<f>{x{n))) = / dfe-"" ^(47rf )- ^ exp 
Jo 



\X{Tn, 



{n)f 



AT 



(14) 
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The exponent was considered part of the worldhne lagrangian for the loop 
path integral, and absorbed into the bosonic worldline Green function. This 
results in a modified Green function 



G, in,r,) = Gs{n,r,) + - f + GsM 

(15) 

valid for Wick contractions of operators inserted into the loop path integral. 

The procedure generalizes to the case of m propagator insertions, and 
leads to modified Green functions 



G'S"'(ti,T2) = GB{n,r2) 

m 

+ 9 [G's(ti,T„J - G's(Ti,TfeJ]At7'[G's(T2,T„J - G's(T2,TfeJ], 

(16) 

with a symmetric m X m - matrix A'™' defined by 



Tki 
Bki 



2 

GB{Ta^,Tai] 



GB{Tak,ni) - Gsin^.Tai) + GB{Tb^,ni] 



;i7) 



(Ti, . . . Tm are the proper-time variables for the inserted propagators). 
4. QED Multiloop Amplitudes 

For scalar electrodynamics, it is obvious what should replace the propaga- 
tor insertion eq.( 13). In this case, the one-loop path-integral eq.( 1) reduces 
to 



t°° di ■> t /'-' /i \ 

V[A] = —e-^'^JVxexp^-J^ dr^-i' + leA^i'')^ (18) 

(note the deletion of the global factor of —2, which takes care of statistics and 
degrees of freedom). This expression may also be interpreted as a Wilson loop 
expectation value. It is well-known, however (see e.g. [18, 20, 21]), that the 
first-order correction to a scalar Wilson loop (due to exchange of one internal 
photon) may be written in terms of a worldline current-current interaction. 
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with A = f - 1. 

For scalar electrodynamics, we will therefore insert one copy of this ex- 
pression into the path integral eq.( 18) for every internal photon. The de- 
nominators will be written in the proper-time representation eq.( 14), which 
leads to the same generalized bosonic two-loop worldline Green function as 
in the scalar case; the numerators will remain, and participate in the Wick 
contractions. 

As in the one-loop case, the transition to spinor electrodynamics may 
then be accomplished by supersymmetrization, which replaces eq.( 19) by 



m. r ,r^,ej' ir^^^l^, . (20) 

Itt^+i Jo Jo ((y - YA^V 



2 47r^+iJo Jo {{X^-Xb) 

In components the double integral reads, after a bit of algebra, 

" '\7f 77 ~ 7f 77^' 



'^ 7 ^-16A(A + 1) . 

[[Xa - Xb) j [[Xa - Xb) j ) 

(21) 

The simplest way to verify the correctness of this naive supersymmetrization 
is to write the one-loop two-photon amplitude in the superformalism, and 
then sewing together the external legs to create an internal photon, using 
Feynman gauge. 

The denominator of eq.( 20) being bosonic, we can again use the proper- 
time representation eq.( 14) to get it into the exponent, and then absorb 
this exponent into the worldline superpropagator. The algebra is completely 
identical to the scalar case, and leads to modified superpropagators 
which are given by the same formulas as in eqs.( 15) and ( 16), with all the 
one-loop Green functions appearing on the right-hand sides replaced by the 
corresponding one-loop superpropagator eq.( 9). The same applies to the 

D_ 

determinant factor (DetA^™')^ , which gives the ratio of the free Gaussian 
path-integral with m propagator insertions compared to the free one-loop 
path integral [14]. 

To summarize, we can obtain a parameter integral representation for the 
sum of all diagrams with one spinor loop and fixed numbers of photons, N ex- 
ternal and m internal, by Wick contracting N vertex operators eq.( 8) with m 
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factors of Jq dTadOa Jq driydOi^D D X^, ^, using the modified superpropagator 
5. The 2-Loop QED /? - function 

As an illustration, we will use this calculus for a rederivation of the two- 
loop QED /3-function, both for scalar and for spinor electrodynamics. We 
will work in component formalism for transparency. 

As usual, matters much simplify if one is only interested in the /3-function 
contribution, as opposed to calculation of the whole amplitude. The simplest 
thing for us to do is to use the effective action formalism with a constant 
background field F^,^, and read off the /3-function from the coefficient of the 
induced F^^F^'^-term. 

As a warm-up, let us first redo the one-loop calculation [5, 7]. 
For constant F^,^, we can choose a gauge such that = ^x^Fpf^. Using this 
A - field in the one-loop path-integral eq.( 1) and expanding the interaction 
exponential to second order, one obtains 



r(i)[F] 



x(-y) £dn dT2^-i>tF,,x\i'iF^px^^ +ViF,:.rii^^F^pi^^^ 
_[47rr]--e-™^^y^ dnj^ dT,[GB\2-Gl,,]j dxoF,^F 



(we have abbreviated Gb{ti, T2) by Gbi2 etc.). The double parameter integral 
(which is really only a single one, as translation invariance can be used to 
set Ti = 0) gives 

£ dn £ dT2 [GB{n,T2f - GF{n,T2f] =-It' . (23) 

To extract the divergent part of the remaining global proper-time integral, 
various regularization methods could be employed. Choosing dimensional 
regularization, we obtain 

/ — e"™ ^ ^ / — e"™ ^T^-T (24) 

Jo T Jo T e 

with e = D — A. Putting things together, the desired one-loop contribution 
to the effective action becomes 

V^'\F]^'^{A7^)-'-^e' j dxoF,^F^^ . (25) 



7 



From this one obtains the one-loop photon wave-function renormahzation 
factor 

(23-l)<" = |-- , (26) 
leading to the usual value for the one-loop QED /? - function, 

= \^ (27) 

O TT 

Now let us describe the two-loop calculation. In the Feynman diagram 
calculation (see e.g. [23]), one would have to separately calculate the three 
diagrams of fig. 1 using some regularization, say dimensional regularization, 
and then extract their ^ - poles. Cancellation of the p- - poles would be 
found in the sum of the results, indicating a cancellation of subdivergences 
due to gauge invariance. 




Figure 1: Diagrams contributing to the two-loop vacuum polarization 

It will be seen that, in the present formalism, the three diagrams get com- 
bined into one calculation. 

Let us begin with the purely bosonic contributions, which correspond to 
the scalar QED calculation. Those are obtained by inserting the worldline 
current-current interaction term eq.( 19) into the bosonic one-loop path- 
integral. After exponentiation of the denominator and absorption into the 
worldline Green function, this results in 

/>oo AT' „ />oo _ rT pT _ D 

= -2(47r)-^/ — e-™^r-T/ df drj dn[f + GsabV' 
Jo 1 Jo Jo Jo 

x(^) I dn I dT2 J dxo-iy'^F.^y'^y^F^fsy^y^ykx) . (28) 
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_D_ 

Note the appearance of the two-loop determinant factor [T + Gs(Ta,Tb)] ^ . 
The Wick contraction of 

kii'xVlvlv^2vlvhx) (29) 

has now to be done, using the two-loop Green function eq.( 15) (care must 
be taken with Wick contractions involving yaiVh-, as the derivatives should 
not act on the Ta,Tb explicitly appearing in that Green function). 

Due to the symmetries of the problem, there are only two nonequivalent 
contraction possibilities. The result is written out in terms of the bosonic 
one-loop Green function and its derivatives. As in the one-loop calculation, 
one next eliminates all factors of Gb appearing by partial integrations with 
respect to Ti,T2,Ta,n. 

As the next step, all fermionic contributions are included by applying the 
one-loop substitution rule, replacing, for example, 

GBl2GB2lGBabGBba {Gb12Gb21 — G FI2G F2l){G BabG Bba — GpabGEba) 

(30) 

etc. 

The integrations must then be carried out. At this stage, what we have is 
the desired contribution to the two-loop effective action in form of a sixfold 
integral (see fig. 2), 

X / dT,dTbdndT2P{T,f,T,,Tb,n,T2)F^,F^'' . 

Jo 

(31) 

The function P (whose bosonic part Pbos is given in appendix A) is a polyno- 
mial in the various Gsij, Gsij, Gpij-, multiplied by powers of [T + Gb {^a, Tfe)] 
In particular, the integrations over ti,T2 are polynomial, and can be per- 
formed easily, either by computer, or using a set of relations of the type 



Jo 



du2Gi2GB23 



2G 



Bl3 



f 

Jo 



du2GBuGB 23 



:G 



B13 + 



1 

30 



(32) 



9 




Figure 2: Definition of the six integration parameters 



which may be derived from the master identities 



Jo 

f 

Jo 



du2 . . . dUnGB12GB2Z ■ ■ - Gb n(n+l) 



du-i . . . du„GF^■2GF■2?. . . .G 



— r5„(|ui - u„+i|)sign"(ui -Un+i] 

2n-l 



n'-J F12'-J F2Z ■ ■ ■ ^Fn(n+1) 



fn- 1)! 



En-i{\ui - u„+i|)sign"(ui - Un+l) . 

(33) 



In writing those identities, we have scaled down to the unit circle again 
(xi = Tui). Bn denotes the n-th Bernoulli-polynomial, and En the n-th 
Euler-polynomial. Due to the fact that those polynomials can be rewritten 
as 



Bn{x) = Pn{x^ — x) (n even) 

B4x) = P4x^-x){x-^) (n odd) 

(34) 

with another set of polynomials Pn{x) (the same property holds true for 
En{x)), the right hand sides can always be reexpressed in terms of Gb,Gb 
and Gp, so that explicit u^'s will never appear in those relations. Those 
integrals needed for the present calculation are listed in appendix B (for 
the spinor-loop case); a proof of the general identities eqs.( 33) is given in 
appendix C ^. 

Next we perform the T - integration, which is trivial: 

■"^For the traces of the left hand sides, recursion relations had already been derived in 

[6]. 
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poo _ _ D I, Ct^ ^ ^ 

/ df[f + GBa,V~~' = ^j±— (A; = 1,2). (35) 
Jo ^ + fc — 1 



Collecting terms, we get 



(•OO _ p 1 p 1 _ 

/ dT dri I dT2 P{T, T, Ta, n, n , T2 
Jo Jo Jo 

-( 



^\(D-^){D-1)g'Jt + {D-2){D-1)gI-2] . 



(36) 



Let us also give the corresponding expression for scalar QED, which is ob- 
tained by using only the bosonic part Pi^os of the function P: 



/ df dn dT2Pko.{T,f,T,,n,n,T2) = -{D-l)Gs:,T 

Jo Jo Jo 6 

32 1 D 1 , D 

HD - 1)(^ - 4)G'^J T + —{D - 2){D - 'J)g'^J 



(37) 



Setting Ta = 0, the integration over Tf, produces a couple of Euler Beta- 
functions, 

dTa I dn Gtal = B{k + l--,k + l--) T'+'-T . (38) 

As in the one-loop case, the remaining electron proper-time integral just gives 
a r - function: 

r ^e-™'^r^-^ = r(4 - i))m2(^-^) . (39) 
JO T 

Combining terms and performing the e - expansions for the effective la- 
grangians, we obtain 



(40) 
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So far this is a calculation of the bare regularized effective action. What 
about renormalization? The counterdiagrams due to electron wave function 
and vertex renormalization need not be taken into account, as they cancel 
by the QED Ward identity [Zi = Z2). However, we have used the electron 
mass as an infrared regulator for the electron proper-time integral eq.( 39); 
mass renormalization must therefore be dealt with. 

Generally, we do not know, at present, how to perform renormalization 
completely in terms of worldline concepts; we have to refer to standard field 
theory for this part of the calculation. 

Our calculation corresponds to a Feynman calculation in dimensional 
regularization and Feynman gauge, so we need to know the corresponding 
one-loop mass renormalization counterterms, both for scalar and spinor QED. 
This is a simple textbook calculation, of which we give the result only: 



scat { A ^ 

— — e (47rj 



seal 



^f^spin 6 2/, \-2 



-e'(47r) 



e 



(41) 



Insertions of those counterterms into the one-loop path integral produce the 
following contributions to the two-loop effective lagrangians. 



' dm 



2 



d 

dm 



AT^liF] = Sm,,,^—Till[F] 



e 



(42) 



Here r^J,'^ denotes the one-loop path integral eq.( 22), F^^^^ its scalar QED 
counterpart. 

Extraction of the f3 - function coefficients proceeds in the usual way. From 
the total effective lagrangians 
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(43) 



one obtains the two-loop photon wave-function renormalization factors, and 
from those the standard results for the two-loop f3 - function coefficients 
[24, 25], 

&H=PiXH = ^ ■ (44) 

Observe that in the spinor-loop case, the integrand after performance 
of the first three integrations, eq.( 36), has only one term which would be 
divergent for D = A when integrated over Tf,. Moreover, the coefficient of 
this term vanishes for D = A. This suggests that this calculation can be 
further simplified by using some four-dimensional regularization scheme. And 
indeed, if we do the spinor-loop calculation in four dimension, then instead 
of eq.( 36) we find simply 

/ df dn dT2P{T,f,T^,n,n,T2) = -8 . (45) 

JO Jo Jo 

This time there is no dependence on Ta,Tb left, so that one immediately gets 

,00 dT 



(2) 
spin 



[F] = {A7r)-'e' j^^ —e-"^'^F,^F^- . (46) 

It is only the final electron proper-time integral that now needs to be 
regularized. This can be done by introducing a proper-time cutoff Tq at the 
lower integration limit, which replaces eq.( 39) by 



°° '^^e-™'^ ~ -Hm'To) (47) 



To T 



(Pauli-Villars regularization could be used as well, though proper-time reg- 
ularization appears more natural in the worldline formalism). With this 
regulator, the two-loop effective lagrangian becomes 

^Zl[F] ~ -Hm'To){A7r)-'e'F,^F^^ + finite . (48) 

In spite of the apparent suppression of subdivergences, there is again a contri- 
bution from mass renormalization, which can be determined by comparison 
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with the corresponding Feynman calculation ^. On-shell renormalization of 
spinor QED using a proper-time cutoff has been studied in refs. [27, 28]. It 
leads to a one-loop mass renormalization counterterm 

— = 31n(m2ro)e2(47r)-' + finite . (49) 
m 

Insertion of this counterterm into the one-loop path integral gives 

om 

~ 21n(m2ro)(47r)"^e^ j dxoF^^F^"" + finite, 

(50) 

so that mass renormalization now just amounts to a sign change for the 
effective lagrangian: 

^Zl[F] + AC%l[F] ~ ln(m^ro)(47r)-^e^F,.F-^ . (51) 
The extraction of the (still scheme-independent) f3 - function coefficient 
Pspini'^) is again standard [28], and leads back to eq.( 44). 

Proper-time regularization could, of course, also be applied to the corre- 
sponding scalar QED calculation. However, here one either has to regulate 
both the electron and the photon proper-time integrals, or to switch to Lan- 
dau gauge, where the Tf, - integral again becomes finite in D = 4 (as we have 
verified). 

A word of explanation may be in place for our use of the one-loop sub- 
stitution rule eq.( 30) in the two-loop context. The continued validity of the 
substitution rule at the multiloop-level is a consequence of the compatibil- 
ity of the superfield formalism with the mentioned sewing-procedure, as will 
be explained in more detail elsewhere [32]. As an explicit check, we have 
performed this calculation in yet another way, namely by writing the one- 
loop four-photon amplitude in the Bern-Kosower representation, and then 
sewing together two of the photon legs to create the two-loop vacuum polar- 
ization amplitude. Albeit requiring considerably more work, this procedure 
ultimately yields exactly the same parameter integrals as the use of the two- 
loop worldline Green function. 

Moreover, the calculation in the dimensional scheme has been checked in 
detail [29] against a Feynman calculation using the second-order Feynman 
rules of [30]. 

^Mass renormalization had been omitted in an earlier, incorrect version of this calcu- 
lation [26]. 
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Let us summarize the properties of this calculation: 

i) Neither momentum integrals nor Dirac traces had to be calculated. 

ii) Only simple one-dimensional integrals were encountered. 

iii) The three diagrams of fig. 1 were combined into one calculation (in 
fact, in this formalism it is somewhat easier to compute the sum than any 
single one of them). 

iv) In the spinor-loop case, introduction of a regulator could be avoided, 
except for the final electron proper-time integration. 

Property iii) is a general property of the formalism, which should become 
increasingly important at higher orders. Property iv) may well be accidental 
to the two-loop far as Feynman gauge is concerned. However, we 

expect it to hold true in general for the recursively determined gauge where 
Zi = Z2 = 1 [31]. 

6. Conclusions 

To conclude, we have shown that the generalization of the Bern-Kosower 
formalism proposed in [14] holds considerable promise as a tool for multiloop 
calculations in quantum electrodynamics. It allows to write down compact 
integral representations combining all Feynman diagrams with one spinor- 
loop and a fixed number of internal and external photons (the extension 
to an arbitrary number of spinor-loops is straightforward, as will be shown 
elsewhere [32]). This kind of sum of diagrams, however, is known to be 
afflicted with extensive cancellations between diagrams. Those cancellations 
are clearly related to gauge invariance, and to the fact that the Feynman 
diagram calculation splits a gauge invariant amplitude into non-gauge in- 
variant pieces. For instance, it is a well-known fact that, at any fixed loop 
order, the higher-order poles cancel and only the ^ - pole persists in the 
sum of all single-spinor loop Feynman diagrams contributing to the vacuum 
polarization. Moreover, individual Feynman diagrams in this sum contribute 
transcendental numbers to the /? - function coefflcients, which happen to can- 
cel out for those few coefflcients which have been calculated [33, 34]. While 
the cancellation of higher-order poles is a well-understood consequence of 
gauge invariance [31], no convincing explanation has been given, so far, for 
the cancellation of transcendentals ^. We hope that the formalism developed 
in this paper will serve to shed new light on this old problem. 
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^Very recently, concepts from knot theory have been used to establish a link between 
both types of cancellations [35]. 
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Appendix A 



We give here the purely bosonic part Pbos of the function P{T, T, r^, Tfe, ti, T2) 
appearing in eq. (31). The full P is obtained from this by application of the 
substitution rule eq.( 11). In writing this polynomial, we have used symmetry 
with regard to interchange of ti and T2 to combine some terms, and omitted 
some terms which are total derivatives with respect to / dri or / dT2 (those 
terms are easy to identify at an early stage of the calculation). 

Pbos = 7^ {^"^iG^ab^BU + ^'^lGBabGB12GBlaGB2b 

+ 8jGBlaGBabGB12[GB2a — GB2b] 

— 4:jGBlaGB2b[GBla — GBlb][GB2a — GB2b] 

+ 187^(7^^;, (7si2[Gsia — GBlb][GB2a — GB2b]^ ■ 

(A.l) 

We have defined 7 = [T + GBab] ^ ■ 



Appendix B 

Integrals occuring in the calculation of the two-loop spinor-QED /3-function: 
I'du, C du,{Gl,, - Gl,,) = -I (B.l) 

JO Jo O 

t dur[G 
Jo 



1 



^313 — Gbu]^ — 2^*^34 (B.2) 



/ dui / du2{GB12GB23GB34:GB41- 

Jq Jo 



82 4 

-G FI2G F23G F3aG FAl) = "2^*534 ~ 2^*534 (B-3) 



^341- 



/ dui / du2{GB13GB32GB2AGi 

Jo Jo 

8 8 

— Gf13Gf32Gf24:Gf4:i) = 4G^34 + -Gb3A — ~ (B.4) 



3 9 
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Appendix C 



To prove the first of the identities in eq. (33), observe that, by construction, 
^Gb is the integral kernel inverting the first derivative ds acting on periodic 
functions. We may therefore write 

Kn{ui — Un+l) := du2 . . .dUnGB12GB2Z- ■ -GBnin+l) 



10 

2" < Ml I ^s"" I U„+i > . 



(C.l) 



This leads to the recursion relation 

^KJu -u') = T' <u \ I u' >= 2Kn-i(u - u'). (C.2) 

du 

We want to show that the same recursion relation is fulfilled by the polyno- 
mial K„ 



K4u - u') := -^B4\u - u'\)sign^{u - u'). (C.3) 



Explicit differentiation yields 



d ~ 2" 

—Kn{u - u) = -B'^{\u - u'\)sign{u - u')sign^{u - u') 

ou n\ 

on 



[n 

2kr,-l{u-u'). 



(C.4) 



Here the recursion relation for the Bernoulli polynomials was used, B'n{x) = 
nBn-i{x). An additional term arising by differentiation of the signum func- 
tion for n odd can be deleted due to the fact that 

S{x)B4\ x\) = S{x)B40) = (C.5) 

for n odd, n > 1. The proof is completed by checking that the master identity 
works for n = 1 [Bi[x) = x — |), and on the diagonal ui = Un+i for any n 
(this special case has already been proven in [7]). 

The proof of the second master identity proceeds in a completely analogous 
way. 



17 



References 

[1] Z. Bern and D. A. Kosower, Phys. Rev. Lett. 66 (1991) 1669; 
Z. Bern and D. A. Kosower, Nucl. Phys. B379 (1992) 451. 

[2] Z. Bern and D. C. Dunbar, Nucl. Phys. B379 (1992) 562. 

[3] Z. Bern, L. Dixon and D. A. Kosower, Phys. Rev. Lett. 70 (1993) 2677 
(hep-ph/9302280). 

[4] Z. Bern, D. C. Dunbar and T. Shimada, Phys. Lett. B312 (1993) 277 
(hep-th/9307001); 

D.C. Dunbar, P. S. Norridge, Nucl. Phys. B433 (1995) 181 (hep- 
th/9408014). 

[5] M. J. Strassler, Nucl. Phys. B385 (1992) 145. 

[6] M. J. Strassler, SLAC-PUB-5978 (1992)(unpubhshed). 

[7] M. G. Schmidt and C. Schubert, Phys. Lett. B318 (1993) 438 (hep- 
th/9309055). 

[8] D. Cangemi, E. D'Hoker and G. Dunne, Phys. Rev. D51 (1995) 2513 
(hep-th/9409113). 

[9] D. Fliegner, M. G. Schmidt and C. Schubert, Z. Phys. C 64 (1994) 111 
(hep-ph/9401221). 

[10] K. Roland, Phys. Lett. B289 (1992) 148; 

G. Cristofano, R. Marotta and K. Roland, Nucl. Phys. B392 (1993) 
345. 

[11] D.G.C. McKeon, Ann. Phys. 224 (1993) 139; 

D.G.C. McKeon, A. Rebhan, Phys. Rev. D48 (1993) 2891. 

[12] C. S. Lam, Phys. Rev. D 48 (1993) 873 (hep-ph/9212296); 
C. S. Lam, Can. J. Phys. 72 (1994) 415 (hep-ph/9308289). 

[13] G. Siopsis, Phys. Lett. 342 (1995) 224 (hep-ph/9410206). 

[14] M. G. Schmidt and C. Schubert, Phys. Lett. B331 (1994) 69 (hep- 
th/9403158). 

[15] E. S. Fradkin, Nucl. Phys. 76 (1966) 588. 

[16] L. Brink, P. Di Vecchia and P. Howe, Nucl. Phys. B118 (1977) 76. 



18 



[17] E. S. Fradkin and A. A. Tseytlin, Phys. Lett. B163 (1985) 123. 
[18] A. M. Polyakov, Gauge Fields and Strings, Harwood 1987. 
[19] J.W. van Holten, Z. Phys. C66 (1995) 303 (hep-th/9408027). 
[20] Yu.M. Makeenko and A. A. Migdal, Nucl. Phys. B 188 (1981) 269. 
[21] H. Dorn, Fortschr. Phys. 34 (1986) 11. 

[22] L. Brink, S. Deser, B. Zumino, P. Di Vecchia and P. Howe, Phys. Lett. 
B64 (1976) 435. 

[23] C. Itzykson and J. Zuber, Quantum Field Theory, McGraw-Hill 1980. 

[24] R. Jost and J.M. Luttinger, Helv. Phys. Acta 23 (1950) 201. 

[25] V. L Ritus, Zh. Eksp. Teor. Fiz. 73 (1977) 807 [Sov. Phys. JETP 46 
(1977) 423]. 

[26] M. G. Schmidt and C. Schubert, HD-THEP-94-32 (hep-ph/9408394), 
and the original version of this preprint. 

[27] W. Y. Tsai, Phys. Rev. D 10 (1974) 1342. 

[28] W. Dittrich and M. Renter, Effective Lagrangians in Quantum Electro- 
dynamics, Springer 1985. 

[29] A. Laser, M. G. Schmidt and C. Schubert, in preparation. 

[30] A. Morgan, Phys. Lett. B351 (1995) 249 (hep-ph/9502230). 

[31] K. Johnson, R. Willey, and M. Baker, Phys. Rev. 163 (1967) 1699. 

[32] M. G. Schmidt and C. Schubert, in preparation. 

[33] J. L. Rosner, Ann. Phys. 44 (1967) 11; 

E. de Rafael and J.L. Rosner, Ann. Phys. 82 (1974) 369. 

[34] S. G. Gorishny, A. L. Kataev and S. A. Larin, Phys. Lett. B 273 (1991) 
141. 

[35] D. Broadhurst, R. Delbourgo and D. Kreimer, OUT-4102-60 (hep- 
ph/9509296). 



19 



